Why do millisecond pulsars have weaker magnetic 
fields compared to ordinary pulsars? 
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Millisecond pulsars, with magnetic fields weaker by three to four orders com- 
pared to those of ordinary pulsars, are presumed to be neutron stars spun up 
by binary accretion. We expect the magnetic field to get screened by the ac- 
creted material. Our simulation of this screening mechanism shows, for the 
first time, that the field decreases by a purely geometric factor sin 7 ^ 2 6 P ;i be- 
fore freezing to an asymptotic value, where 6^,i is the initial angular width of 
the polar cap. If Pji lies in the range 5°-10°, then the magnetic field diminution 
factor turns out to be ~ 10 3 -10 4 in conformity with observational data. 
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Pulsars — which are rotating neutron stars — are seats of strongest magnetic fields known to 
mankind. Magnetic fields of most pulsars are around 10 n -10 13 G, whereas the typical rotation 
periods are about 1 s. However, there exist a handful of known pulsars with considerably smaller 
rotation periods, which also have much weaker magnetic fields around 10 8 G. These so-called 
millisecond pulsars are often found in binary systems (0). 

It is believed that most of the neutron stars (barring possible magnetars ©) are born with 
characteristics typical of ordinary pulsars. If the neutron star happens to be in a binary system, 
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then it is possible for it to accrete matter with angular momentum from the binary companion. 
Millisecond pulsars are thought to be neutron stars which have been spun up in such a binary 
accretion process (3). Since their magnetic fields are weaker by a factor 10 3 -10 4 compared 
to the magnetic fields of ordinary pulsars, presumably the magnetic field of the neutron star 
also decreases during the accretion phase. Several alternative scenarios have been proposed to 
explain how this decrease of magnetic field takes place, starting from models relating the field 
evolution to the spin evolution to models based on accretion-induced Ohmic decay in the heated 
crustal layers (@|5||6][7|0l2llZIl- One possibility is that the magnetic field gets buried under the 
accreted matter. Although this idea has been around for a while (ETHT^HZHIZHIZSllZSllZZh we 
have, for the first time, carried out a detailed 2-D simulation to test this idea. Details of our 
calculations are presented elsewhere (|7S]|Z3- One of the attractive features of this scenario is 
that the factor 10 3 -10 4 by which the magnetic field decreases can be explained in a very elegant 
and simple way as arising purely out of geometric considerations. 

The strong magnetic field of an accreting neutron star channelizes the accreting material 
to flow through the polar regions. As the magnetic field of the neutron star decreases because 
of the screening due to accreting material, it is less able to channelize the accretion flow and 
thereby the polar cap widens. One can easily find out how the angular width 6 P of the polar 
cap depends on the surface magnetic field B s of the neutron star (see, for example, Shapiro and 
Teukolsky (12771)). The field line starting from 9 P at the surface of the neutron star, with a radius 
r s , is the last closed field line of the dipolar field and passes through the Alfven radius r A . It 
easily follows that 



Assuming that the ram pressure of the freely in-falling accreting material at the Alfven radius 
equals the magnetic pressure, a few steps of easy algebra give 




(1) 



r A = {2GM)- 1 ' 7 rl 2 ' 7 BfM 
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(2) 
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where M is the mass of the neutron star and M the accretion rate. It follows from (1) and (2) 
that 

sintfp oc B; 2/7 . (3) 

This is how the polar cap widens with the weakening magnetic field until 9 P becomes equal to 
90° when (3) obviously ceases to hold. On taking M = 10 33 gm, M = 1(T 8 M yr _1 , r s = 10 
km, B s = 10 12 G, we find from (2) that r A ~ 300 km. Substituting this in (1), we conclude that 
the initial polar cap angle is of order 10°. 

The accreting materials falling through the two polar caps flow horizontally towards the 
equator in both the hemispheres. At the equator, the opposing materials flowing in from the 
two poles meet, sink underneath the surface (inducing a counter-flow underneath the equator- 
ward flow at the surface) and eventually settle radially on the neutron star core. With a suitably 
specified flow having these characteristics, we have studied kinematically how the magnetic 
field evolves with time, taking into account the fact that the polar cap width changes with the 
evolution of the magnetic field, thereby altering the velocity field also. Fig. 1 shows the velocity 
field at an early stage (A) and at a late stage (B). We find that the equator- ward flow near the 
surface seen in Fig. 1A is quite efficient in burying the magnetic field underneath the surface. 
However, when the polar cap opens to 90°, the accretion becomes spherical and radial, as seen 
in Fig. IB. It is found that such accretion is not efficient in burying the magnetic field any 
further. Fig. 2 shows magnetic field lines at different stages of evolution, whereas Fig. 3 plots 
the surface magnetic field at 45° as a function of time. Clearly, the magnetic field at the surface 
of the neutron star keeps decreasing until the polar cap opens to 90°, after which the magnetic 
field is essentially frozen, since the radial accretion cannot screen it any further. If 9 Pii is the 
initial polar cap width, then it follows from (3) that the magnetic field would decrease by a 
factor (sin 90°/ sm9 P - 1 ) 7 ^ 2 from its initial value before it is frozen to an asymptotic value. On 
taking 9 P :i in the range 5°-10°, this factor turns out to be about 10 3 -10 4 , exactly the factor by 
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which the magnetic fields of millisecond pulsars are weaker compared to the magnetic fields of 
ordinary pulsars. 

Put another way, the magnetic field freezes when the Alfven radius becomes equal to the 
neutron star radius. The asymptotic value of the surface magnetic field can be found directly 
from (2) by setting ta equal to r s , which gives 

5 asymp = (2GM)^M 1/2 r^\ (4) 

On using the various standard values mentioned before, we find -B a sym P ~ 10 8 G. When the 
magnetic field falls to this value, it can no longer channelize the accretion flow, resulting in the 
flow becoming isotropic. Such a flow is unable to screen the magnetic field any further. After 
the accretion phase is over, the neutron star appears as a millisecond pulsar with this magnetic 
field. We propose this as the reason why millisecond pulsars are found with magnetic fields of 
order 10 8 G. 

We would like to thank Somnath Bharadwaj and Sayan Kar for their helpful comments. 
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Fig. 1. Velocity fields induced inside the neutron star due to accreting material settling on the 
surface. The arrows indicate the flow velocities, whereas the dashed lines are contours of con- 
stant V. (pv) . The velocity field is specified in the code in such a way that it keeps changing with 
the evolution of the magnetic field. The details of how we do this are provided elsewhere (I7P1) . 
Here we show the velocity fields at two different instants: (A) when the magnetic field is strong 
and the polar cap angle is small (at a relatively early stage); and (B) when the magnetic field has 
become much weaker and the polar cap has opened up (at a relatively late stage). The dashed 
lines (contours of constant V.(pv) ) indicate regions which are sources of new material due to 
accretion. We see in (A) that the new material is dumped in a narrow polar cap, inducing an 
equator-ward flow just below the surface. After reaching the equator where this flow meets the 
oppositely-directed flow from the other pole, the flow sinks underneath the surface, induces a 
counter-flow and eventually settles on the core of the neutron star. On the other hand, we see in 
(B) that the new material falls isotropically all over the surface and the induced flow is radially 
inward. The horizontal flow shown in (A) is expected in reality to be confined only within 1% 
of the neutron star radius immediately below its surface (18). Here we have inflated that layer 
to 10% of radius for easy visualization. 

Fig. 2. Magnetic field lines during different phases of evolution. The dashed lines indicate the 
contours of constant V.(pv) at the same instants. The topmost panel shows the initial magnetic 
field, which is evolved by solving the induction equation numerically, with the specified veloc- 
ity field which keeps changing as the magnetic field weakens. (The velocity fields at an early 
stage and at a late stage are shown in Fig. 1.) The details of the numerical code are given in the 
Appendix of our earlier paper (f7%t. which also discusses the various boundary conditions used. 
The code has been adopted from a code which was used extensively for studying the evolution 
of the solar magnetic fields (I2ill22l) . 
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Fig. 3. We show how the magnetic field at the surface at 45° latitude changes with time (solid 
line). The evolution of the polar cap angle 9 P with time is also indicated (dashed line). It is 
clear the magnetic field decreases rapidly till the polar cap opens to about 90°, after which the 
decay of the magnetic field is significantly halted. 
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Figure 1: 
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Figure 2: 
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